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We use the LoopSim and VBFNLO packages to investigate a merged sample of partonic events 
that is accurate at NLO in QCD simultaneously for the WZ and WZ+jet production processes. 
In certain regions of phase space such a procedure is expected to account for the dominant part 
of the NNLO QCD corrections to the WZ production process. For a number of commonly used 
experimental observables, we find that these corrections are substantial, in the 30 — 100% range and 
in some cases their inclusion can reduce scale uncertainties by a factor of two. As in the underlying 
VBFNLO calculations, we include the leptonic decays of the vector bosons and all off-shell and 
finite-width effects. 
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I. INTRODUCTION 

The study of di-boson production processes at the LHC 
is important both to test the Standard Model (SM) and 
because they constitute relevant backgrounds to the be- 
yond standard physics (BSM) searches. They are for 
example sensitive to tri-linear gauge couplings (TGC) 
fixed by the underlying electroweak gauge group in the 
SM. Any deviation in the values observed at experiments 
would indicate new physics effects. In particular, WZ 
production measurements [l|, 0] are relevant to constrain 
the WWZ tri-linear coupling and to search for charged 
heavy bosons from BSM sector via Jacobian peaks 

To match the experimental accuracy, precise and reli- 
able predictions beyond the leading order approximation 
are required. The next-to-leading order (NLO) QCD 
corrections for WZ production were computed in 0, H[ 
and turned out to be sizable. These large corrections 
are a well known fact for colorless production channels. 
At NLO QCD, new topologies and new partonic sub- 
processes appear, Fig.[TJ e.g. those with a W or Z emitted 
from a jet or an initial-state parton. These new config- 
urations can have partonic luminosities (e.g. qg) that 
are enhanced over those present at LO (e.g. qq), which 
can lead to large corrections in differential distributions. 
They also permit soft or collinear bosons emissions from 
a quark which results in a s a| w In 2 (pTj/my) enhance- 
ment for a number of observables. 

At the next-to-next-to-leading order (NNLO) further 
new sub-processes and topologies, including those with 
soft and collinear bosons, start to contribute, Fig. [TJ 
which might also lead to large corrections in some re- 
gions of phase space. 1 

Therefore, using only the NLO results to compute the 
WZ backgrounds may lead to a misinterpretation of a 
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FIG. 1: Example diagram contributing to WZ production at 
LO, NLO and NNLO. 



possible excess observed in data, which could be wrongly 
attributed to new physics. Thus, it is of great interest 
to try to assess the genuine NNLO QCD contributions 
for this process. The current state of the art is that the 
virtual two-loop corrections for WZ production are un- 
known. However, recently, the NLO QCD corrections 
for WZ+jet were computed in [Ioj |. including anomalous 
couplings [III, and they are available in the VBFNLO 
package [12(- This result provides the mixed real-virtual 
and the double real C(a 2 ) contributions to the WZ pro- 
duction including the NLO corrections to the qg chan- 
nel. This opens the possibility to merge the results for 
WZ@NL0 2 and WZj@NLO and obtain the part of the 
NNLO result for WZ production that is associated with 
channels absent at LO, which is the dominant contribu- 
tion for a range of differential distributions at high px- 

This can be done in a consistent way by employing 
the LoopSim method (l3j which uses unitarity to cancel 
infrared and collinear divergences that arise when one 
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adds tree level results with different multiplicities. The 
LoopSim prescription is general and it has been found 
to work well for Z+jet and Drell-Yan production, es- 
pecially for observables which receive significant correc- 
tions at NLO [13(. The method supplements a set of 
real and real- virtual diagrams with approximate contri- 
butions from loop diagrams. The result at a given order 
is denoted with h to distinguish it from the exact result. 

In this paper, we compute the nNLO corrections (ac- 
cording to the notation introduced above) to WZ pro- 
duction at the LHC 

PP ^ w ± z + x^-if tii+e^ + x, (i) 

including the leptonic decays and full off-shell and finite 
width effects. 

We organize this paper as follows: In SectionQll we give 
technical details of our computation including a short 
description of the LoopSim method. In Section IIII1 we 
present numerical results and the impact of the nNLO 
QCD corrections on various differential distributions. In 
Section [IV! we summarize our findings. 

II. CALCULATIONAL DETAILS 

For the calculation of WZ@nNLO, we used the im- 
plementations of WZONLO and WZj@NLO from the 
VBFNLO package [12], with the latter providing the 
double-real and the real-virtual parts of the NNLO re- 
sult. Then the missing two- loop contributions are ap- 
proximated by LoopSim from the tree-level and one-loop 
parts of WZj@NLO. To allow for communication between 
VBFNLO and LoopSim, we developed an interface be- 
tween the two programs. On one hand, it consists of an 
extension of VBFNLO which enables it to write events 
in the Les Houches (LH) format at NLO and, on the 
other hand, of a class which reads the Les Houches events 
and passes them to LoopSim. Then LoopSim assigns 
an approximate angular-ordered branching structure to 
each LH event, with the help of the Cambridge/ Aaachen 
(C/A) [3, [H| jet algorithm. By default, it uses a ra- 
dius i?LS — 1 an d by varying i?LS; onc probes one class 
of uncertainty of the method. As we shall see in the 
next section, except for the very low pt region, this un- 
certainty is significantly smaller than that coming from 
renormalization and factorization scale variation. 

The error that we make for an observable (A) is a finite 
constant associated with the LO topology and it is free 
of infrared and collinear divergences since these are, by 
construction, suitably treated by LoopSim. This implies 

4n ) lo-4n ) lo=0(«Mo ) )- ( 2 ) 

Thus, differential distributions sensitive to new channels 
and new kinematically enhanced configurations resulting 
in large NLO K-factors, should be close to the full NNLO 
result, providing predictions more precise than the pure 
NLO corrections to WZ. 




H T [GeV] 

FIG. 2: Differential cross sections and K factors for the 
effective mass observable, defined in Eq. Q, for the LHC at 
\fs — 8TeV. The bands correspond to varying /if = fin 
by factors 1/2 and 2 around the central value from Eq. {3J. 
The cyan solid bands give the uncertainty related to the _Rls 
parameter varied between 0.5 and 1.5. The distribution is a 
sum of contributions from two unlike flavor decay channels, 
eefiVfj. and /i/ie^ e . 



III. NUMERICAL RESULTS 

In our computations, regardless of the order, we use the 
MSTW NNLO 2008 (3 PDFs with a s (M z ) = 0.11707. 
We choose m z = 91.1876 GeV, m w = 80.398 GeV and 
Gf = 1.16637 x 10~ 5 GeV~ 2 as electroweak input pa- 
rameters and derive the weak mixing angle from the 
Standard Model tree level relations. The center-of-mass 
energy is fixed to y/s = 8TcV if not specified other- 
wise. As in the underlying VBFNLO WZ@NLO and 
WZjQNLO calculations, only the channels with W and 
Z/7* decaying into unlike flavor leptons are considered, 

i.e. pp -> W ± Z + X -> if z7i + x - From thc 

theoretical perspective, to take into account all possible 
decay channels, i.e. eeev e , ee\xv^ nnev e , fififiv^, one can 
safely multiply the single flavor result by a factor four, 
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FIG. 3: Differential cross sections and K factors for the pt of the hardest lepton for the LHC at ^/s = 8 TeV (left) and 
\fs = 14 TeV (right). The bands correspond to varying /if = by factors 1/2 and 2 around the central value from Eq. (f3|. 
The cyan solid bands give the uncertainty related to the Rls parameter varied between 0.5 and 1.5. The distribution are sums 
of contributions from two unlike flavor decay channels, ee^iv^ and fifiei/ e . 



since the Fermi interferences due to identic al p articles 
in the final state are below per mille level [lfj. How- 
ever, as explained later, we shall only apply a cut on the 
reconstructed mass of the opposite charge same flavor 
lepton pair. Therefore, all the results shown in Figs. [5]- 
[T] correspond to the sum of contributions from the two 
decay channels, ee^iv^ and [i[iev e . All off-shell effects 
are included, which takes into account spin correlations 
and also virtual photons. Finite width effects, related 
to leptonic decays of the electroweak gauge bosons, are 
accounted for with a fixed- width scheme [lOj . Top quark 
effects are not considered and all other quarks are taken 
massless. Effects from generation mixing are neglected 
since the CKM matrix is set to the identity matrix. As 
the central value for the factorization and renormaliza- 
tion scales we choose 



f*F,R = y^2pT,pm-tons+2\JPT,W + m W+2\J P T,Z + m Z 

(3) 



where Pt,v and my have to be understood as the recon- 
structed transverse momenta and invariant masses of the 
decaying bosons. To study the impact of the QCD cor- 
rections, we choose inclusive cuts. The charged leptons 
are required to be hard and central: pr,i > 15(20), for 
I coming from Z(W), and \yi\ < 2.5. The missing trans- 
verse energy must satisfy the cut -Er.miss > 30 GeV. The 
reconstructed mass from the opposite-charge same-flavor 
leptons has to lie in the window 60 < < 120 GeV, 

which also avoids singularities coming from off-shell pho- 
tons, 7* — > . We cluster all final state partons with 
\Up\ < 5 to jets with the anti-fc t algorithm [171 ]. as imple- 
mented in Fast Jet 0,1111, witn tne radius R = 0.45. For 
observables that involve jets, we consider only those jets 
that lie in the rapidity range |j/j e t| < 4.5 and have trans- 
verse momenta prjet > 30 GeV. In addition, we impose 
a requirement on the lepton-lepton and lepton-jet separa- 
tion in the azimuthal angle-rapidity plane ARmj-\ > 0.3. 

As a first check of our setup, we have merged WZ@LO 
and WZjOLO to produce WZ@nLO, which can be tested 
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FIG. 4: Differential cross sections and K factors for the missing transverse energy (left) and the transverse mass of the WZ 
system (right) for the LHC at ^fs — 8 TeV. The bands correspond to varying /if — by factors 1/2 and 2 around the central 
value from Eq. ([3]). The cyan solid bands give the uncertainty related to the -Rls parameter varied between 0.5 and 1.5. The 
distribution are sums of contributions from two unlike flavor decay channels, eefif^ and fj,fj,ev e . 



against the full WZQNLO result. In Fig. [2] we show the 
effective mass defined by 

H T = y^Prjets + y^Pr,; + E T>miss , (4) 

which often enters in super-symmetry searches [2(], [2l[ . 
This distribution is very sensitive to the enhancements 
from soft or collinear emissions of the electroweak bosons 
as well as to additional parton radiation coming from new 
channels. 

In the middle panel of Fig. [21 we show the K factors 
with respect to the LO result and, in the bottom, the 
corresponding ratios to NLO. At low Ht values, the dif- 
ference between nLO and NLO is at the level of 30—40%. 
This is due to the fact that we do not have control over 
the finite virtual constant terms which are missing in our 
fiLO approximation. However, as the Ht increases, the 
fiLO result converges to the full NLO, providing a predic- 
tion more accurate than a simple LO calculation. Note 
in the middle panel the large fiLO corrections and fast 



convergence to the NLO result, yielding K- factors of or- 
der of 10. The nNLO corrections can be as large as 100% 
compared to NLO (bottom panel of Fig. [5]) and they are 
clearly beyond the NLO scale uncertainties. We observe 
that the Rls uncertainties are small and there is only a 
marginal reduction in the scale uncertainties at fiNLO. 
The latter is due to the fact that we are favoring regions 
of the phase space associated with new topologies enter- 
ing at NNLO which are only computed at LO. 

The integrated cross section for one lepton flavor, 
cr(pp — > hv^hh), dominated by the low Ht region, in- 
creases only by about 5% from 25.7 ± 1.0 (scale) fb at 
NLO to 26.9 ± 0.9 (scale) ± 0.4 (R LS ) fb at fiNLO. We 
have also computed the total cross sections for 7 TeV at 
NLO and the result agrees with this quoted in Q. 

In Fig. [3] we show differential distributions for the lep- 
ton with maximum pt for y/s = 8 and ^fs = 14 TeV, 
left and right panels respectively. We observe that the 
fiNLO corrections are large at both energies. For the 
14 TeV run, these corrections are beyond the rcnormal- 
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ization and factorization scale uncertainties above pr = 
150 GeV. The same is true for the 8 TeV case only at 
slightly higher px- Even though the effect is less pro- 
nounced than at 14 TeV, due to the smaller relative im- 
portance of additional parton radiation, the corrections 
are at the 15% level already at 200 GeV. In addition, for 
the 14 TeV case, we show the curves with a veto on the 
jets at NLO and fiNLO where we require the absence of 
any jets with pt > 50 GeV. The nNLO-veto corrections 
are negative and beyond the scale uncertainties of the 
NLO-vetoed predictions. This might be of relevance for 
anomalous coupling searches since vetoed distributions 
are often used to suppress additional radiation. The en- 
tire nNLO-veto result has larger scale uncertainties than 
the NLO-vetoed curves, thus revealing, partially, acci- 
dental cancellation happening at NLO. Such a feature 
has been extensively discussed for jet vetoes also in the 
context of NNLO calculations of Higgs-boson production, 
e.g. [HI and references therein. 

In Fig. 2] (left), we show the differential distributions 
for the missing transverse energy. We see that the fiNLO 
corrections to this observable can be as large as the 30%, 
exceeding the scale uncertainty of the NLO result. In 
the right-hand plot of Fig. [4j we show the cluster energy 
(transverse mass of the WZ system) defined by 

mr.wz = yf(E™ + E^y - (pf +pf )2 - (pW +P zy , 

(5) 

where E T ' and p^ y z are the transverse energy and the 
transverse momentum components of the bosons recon- 
structed from the four-momenta of the leptons. The 
fiNLO corrections to this observable are small, which 
means that this distribution is not particularly sensi- 
tive to the new topologies that appear at NNLO. For 
the same reason, the finite terms from the two-loop di- 
agrams, which are missing in the fiNLO result, are of 
larger relative importance for this observable than they 
are for H?, Pt.i, max or Et, miss: which is also reflected in 
larger i?LS uncertainties in the mT,wz case. 

The reason behind the marked difference in the rela- 
tive size of the fiNLO corrections between m^wz and 
all the other observables which we have discussed is the 
following. For observables like PT,z,max or -Er.miss, to 



receive large contributions in the high-p^ region, it is 
enough that only one of the two bosons is produced at 
high pt- This high-py boson recoils against a high-p^ 
QCD parton. The other boson can be in particular soft 
and collinear to a quark and such configurations, one of 
which is depicted in the bottom diagram of Fig. [TJ come 
with the a^ag W ln 2 (pTj/my) enhancement. In the case 
of tot.wz, however, the favored configurations are those 
in which both bosons have sizable p^s and are prefer- 
ably back to back therefore do not lead to logarithmic 
enhancements at NNLO. 

IV. CONCLUSIONS 

In this letter, we have used LoopSim together with 
VBFNLO to compute an approximation to the NNLO 
QCD corrections for the process pp — > ifpil^i^ + X- 
Our result, referred to as fiNLO, is expected to be ac- 
curate in the regions of phase space dominated by the 
topologies other than those present at LO. As in the un- 
derlying VBFNLO WZONLO and WZj@NLO calcula- 
tions, our treatment includes the leptonic decays of the 
vector bosons and all off-shell and finite-width effects. 

We found that the fiNLO corrections to a number of 
observables are sizable at high px and have non-trivial 
kinematic dependence. It is therefore important to take 
them into account in searches for physics beyond the SM 
and other physics analyses that involve WZ production. 
The VBFNLO+LoopSim code is available from the au- 
thors on request and it will be made public in the near 
future. 
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